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ABSTRACT 

We report the results from spectroscopic observations of 1 13 ultra-wide, low-mass binary systems, composed 
largely of M0-M3 dwarfs, from the SLoWPoKES catalog of common proper motion pairs identified in the 
Sloan Digital Sky Survey. Radial velocities of each binary member were used to confirm that they are co- 
moving and, consequently, to further validate the high fidelity of the SLoWPoKES catalog. Ten stars appear 
to be spectroscopic binaries based on broad or split spectral features, supporting previous findings that wide 
binaries are likely to be hierarchical systems. We measured the Ha equivalent width of the stars in our sample 
and found that components of 8 1 % of the observed pairs has similar Ha levels. The difference in Ha equivalent 
width amongst components with similar masses was smaller than the range of Ha variability for individual 
objects. We confirm that the Lepine et al. ^-index traces iso-metallicity loci for most of our sample of M 
dwarfs. However, we find a small systematic bias in C, especially in the early-type M dwarfs. We use our 
sample to recalibrate the definition of While representing a small change in the definition, the new ^ is a 
significantly better predictor of iso-metallicity for the higher mass M dwarfs. 

Subject headings: binaries: visual — binaries: spectroscopic — stars: abundances — stars: kinematics — 
stars: low mass, brown dwarfs — stars: magnetic fields stars: subdwarfs 



1. INTRODUCTION 

Low-mass stars, generally defined as the regime brack- 
eted by the hydrogen-burning limit (^0.08 M©) and the on- 
set of molecular lines in the photosp here (^0.8 MfT,), make 
up 70% of the Milky Way's stars (iBochanski et al.1 120101) 
and are, perhaps, the best tracers of the structure, dynam- 
ics, and evolutionary history of the Galaxy. However, their 
intrinsic faintness has historically limited the construction of 
large samples. In addition, the ubiquitous molecular features 
in their photospheres and the resulting incomplete line lists 
has restricted the accuracy and usefulness of theoretical at- 
mospheric models. Larg e surveys, such as the Sloan Dig- 
ital Sky Survey (SDSS; Yorket al.''2000^ and the Two Mi- 
cron All Sky Survey (2MASS; Skrutskie et al. 2006), have 
played a large role in advancing our understanding of low- 
mass stars. With a photometric catalog of more than 33 mil- 
lion (jBochanski et al. 2010) and a spectroscopic catalog of 
more than 70,000 (.West et al.. .2011.) M dwarfs SDSS has 
enabled studies of the spatial (IBochanski et al.l 1201 Oh and 
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kinem atic distributions (IBochanski et al.ll2b07bt iFuchs et akl 
|200l in the Milky Way; the mass and lumin osity func- 
tions (iCovev et al.l l2008t IBoch anski et all 120101): and mag- 



netic activity (e.g..lWest et al.ll2 008. 201 iF lKruse e"taLll25Tol: 
iKowalski etaLll2009l; iHilton et al.il2010i) of low-mass stars. 

The metallicity of low-mass stars remains an elusive pa- 
rameter to measure. Given the large number of M dwarfs 
in the Milky Way, an absolute metallicity scale tied to an 
easily observable spectral index would allow for the tracing 
of the format ion history and the chemical evolution of the 
Galaxy (e.g.. IWest et al.l 12008). the dependence of the fun- 
damental mass-radius relati on on metallicity at the bottom 
of the main sequence (e.g., iLopez-M orales 2007), and the 
relations hip between meta l licity and the presence of plan- 
ets (e.g., |L^iitin|200l |wSE&FiichiBI^^ While 
spectral modeling has allowed for metallicity determina- 
tions and well-defined metallicity indices for warmer stars, 
such efforts in the late- K and M spectral type regimes (e.g., 
IHauschildt. Allard. & Ba'ronI [19991: iWitte et [20 ill) have 
met with notable problems due to the onset of broad molec- 
ular lines at <4300 K and due to incomplete molecular line 
lists. Some authors have tried to use photometric indices to in- 
fer the metalUci ty (Bonfils et al. 2005; Johnson & Apps 2009] 
Schlaufman & Laughlin 2010), but these techniques rely on 
trigonometric parallax measurements which are uncommon 
forM dwarfs. 

Some useful spectral features that correl ate with metallicity 
have been identified. In the near-infrared, iRoias-Avala et al.l 
(l2010li) developed a metallicity indicator based on the strength 
of the Na I doublet, the Ca I triplet, and a temperature- 
sensitive water index. This technique has so far only been 
calibrated over a limited range but delivers the greatest 
precision among current techniques. Meanwhile, much 
effort has gone into optical spectra. As the TiO band 
in the optical spectrum becomes weaker with decreas- 
ing metallicity (iBesselll 119821) . the ratio of CaH and TiO 
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molecular bands has been used to distinguish M dwarfs 
from M subdwarfs dKirkpat rick. Henrv. & McCarthy 
19911 iRe id, Hawley. & Gizis' 199 5t: jGiz is 1997: 
Lepine. Ric h. & Shara 2003; Bur gasser & Kirkpatric ra 
2006"). Buildin g on th ese studies, Lepine, Ri ch. & Sharal 
P007, hereafter l LRSp7[) defined the metallicity-dependent 
quantity C using the iReid et al.1 (|i995) CaH2, CaH3, and 
Ti05 molecular band heads; this allowed for the segrega- 
tion of low-mass dwarfs into four classes: dwarfs (dMs), 
subdwarfs (sdMs), extreme subdwarfs (esdMs), and ultra 
subdwarfs (usdMs). These classes may also trace the Galactic 
populations to which these stars belong: dMs were formed in 
the thi n disk, s dMs in the thick disk, and esdMs/usdMs in the 
halo. ILRS07I calibrated the definition of ( using the visual 
binary pairs known at the time , including four sdM and two 
esdM pairs. IWoolf. Lepine. & Wallerstein (2009) mapped the 
( index to an absolute metallicity scale using dM binaries 
with a FGK companion of measurable absolute metallicity; 
but it suffers from significant scatter (^0.3 dex). 

Wide binary (or multiple) systems are ideal, coeval lab- 
oratories to constrain and calibrate the observable proper- 
ties of stars as the components were presumably formed at 
the same time and from the same primordial material but 
have ev olved independently. In lDhital et al.l (2010. hereafter 
iPaper t ) we identified the Sloan Low-mass Wide Pairs of 
Kinematically Equivalent Stars (SLoWPoKES) catalog con- 
sisting of 1342 ultra-wide, low-mass common proper motion 
(CPM) binary system s from the SDSS Data Release 7 (DR7; 
lAbazajian et al.l l2009) by matching angular separations, pho- 
tometric distances, and proper motions. The binary systems 
in the catalog have at least one low-mass (spectral subtype K5 
or later) component, projected physical separations of ^lO'^- 
10^ AU, and distances of -50-800 pc. While most SLoW- 
PoKES pairs are disk dwarfs, 70 low-metallicity sdM and 21 
white dwarf-dM pairs were identified based on their reduced 
proper motions. A Galactic model — ba sed on empirical con- 
straints on the stellar number density (Bochanski et al.ll2()Tol: 
lJuric et al.lf2 008) and velocity dBochanski et al.ll2007al) distri- 
butions in the Milky Way — was used to assess the probability 
that the candidates were a chance alignment of random stars; 
only pairs with such probabilities < 5% were published in the 
SLoWPoKES catalog. The overall fidelity of the catalog is 
expected to be —98%. Hence, the SLoWPoKES catalog is a 
very clean and diverse source of CPM binary systems to be 
used in follow-up studies. As the SLoWPoKES catalog spans 
a wide range in mass and a smaller, but still considerable, 
range in metallicity, it is an ideal sample to constrain the ( in- 
dex as well as to eventually map it to an absolute metallicity 
space. 

Magnetic activity has been shown to decline with age, with 
activity lifetimes of —1-2 Gyr for M0-M3 and —7-8 Gyr for 
M5-M7 dwarfs (IWest et al. 2006, 2008, 2011). This mono- 
tonic decline of activity with age is a signature of stellar spin- 
down and suggestive of a gyrochr onology-like age-rqtation- 
activity relationship in M dwarfs dSkumanichll 19721; iBarnesI 
I2003L I2OO7I; iDelorme et al.ll20l Oi). Leveraging the coevality 
of components of the SLoWPoKES pairs is a good way of 
testing this relationship. 

We have caiTied out a spectroscopic follow-up study of 113 
CPM pairs from the SLoWPoKES catalog. Section |2] details 
our observations and the data reduction procedures. In Sec- 
tion |3] we use our radial velocities to assess the fidelit y of the 
observed SLoWPoKES pairs, use them to redefine the lLRSOTi 

index, and examine the magnetic activity properties of the 



SLoWPoKES pairs. The conclusions are presented in Sec- 
tion S] 

2. OBSERVATIONS & DATA REDUCTION 

The spectroscopic targets were selected from the SLoW- 
PoKES catalog based on their brightnesses, colors, and in- 
ferred mass ratios. Both components were required to be 
brighter than r — 17 so as to obtain the desired S/N within 
a reasonable integration time. Efforts were made to obtain 
(i) an even distribution in r — z space for both the primary 
and secondary components and (ii) a roughly equal number 
of equal-mass (within 5% of each other) and unequal-mass 
(> 5% of each other) systems. We estimated masses from the 
r — z colors based on lKraus & Hillenbrand (2007). 

Observations were carried out with the GoldCam spectro- 
graph on the KPNO 2.1m telescope on two separate observ- 
ing runs on January 1 1-16, 2009 UT and March 26-31, 2010 
UT. For both runs, the #36 grating (1200 lines mm^^) in the 
first order, blazed at 7500 A, along with the OG 550 order- 
blocking filter were used resulting in a wavelength coverage 
of -6200-8200 A with a dispersion of 0.62 A pixel" 1. A slit 
width of 2" was used to maximize the number of photons col- 
lected yielding an effective resolution of 1 .8 A and a resolving 
power of 3500. Both components of a binary were observed 
at the same time by rotating the slit to align with the position 
angle of the binary. While the rotation had to be done man- 
ually and required —10 min of overhead time, it was more 
efficient than observing each component separately. 

Each night quartz flats and biases were taken before the tar- 
gets were observed; when the first half of the night was lost 
due to weather, the flats and biases were taken in the morning. 
For wavelength calibration, HeNeAr comparison arcs, along 
with the BG 38 order-blocking filter, were generally taken af- 
ter each target or when the CCD orie ntation was rota t ed. A 
suite of radial velocity standards from iDelfosse et al] (119981) 
were observed, which we used to assess our radial velocity 
precision (see below). Similarly, a flux standard — HZ 44, a 
bright sdO star — was observed each night during the second 
run. Both observing runs were conducted in bright time, of- 
ten during non-photometric seeing. A combination of clouds 
and high winds caused the loss of 3-4- nights between the two 
runs. 




-100 -50 50 100 
(Delfosse et al.) 

Figure 1. Comparison of the radial velocities for the standard stars 
from the 2009 (triang les) and 2010 (pluses) observation runs selected 
from lDelfosse et aiTfT998i) shows no system atic trends. The m.a.d. 
of the difference between the lDelfosse et all {1998) and the measured 
Vr was 5.7 km s~^; we adopt this value as the error in our radial 
velocity measurements. The 3 m.a.d. regime is shaded in gray. 
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Table 1 

Radi al Velocity Standards from 
jPelfosse et aU1998l) 



Name Spectral Type Vr (km s ^ ) 



GJ 1057 


M4 


27 


GJ 1093a 


M4 


-30 


GJ 1 1 1 1 


M8 


9 


GJ 1156 


M5 


4 


G170 


M2 


-26 


Gl 105b 


M3 


26 


Gl 109 


M2 


30 


G 165-08 


M4 


8 


G1205 


MO 


8 


G1251 


M2 


22 


G1338 


MO 


11 


Gl 380 


K5 


-26 


G1411 


M2 


-85 


G1412B 


M5 


68 


G1450 


Ml 





G1514 


MO 


14 


GI581 


M2 


-10 


G1623 


M2 


-27 


G1625 


Ml 


-13 


LHS 1805 


M4 


1 


LHS 2520 


M3 


80 


LHS 1885 


M4 


16 



All spectra were bias-subtracted, flat-fielded, extracted, co- 
added, wavelength-calibrated, and flux-calibratecQ using stan- 
da rd IRAlf^ p ro cedure s, following the prescription detailed 
in iMassev et alj ( 119921) . Eighteen pairs where one of the 
components (usually the fainter secondary) was not well- 
calibrated or had low S/N were removed from the sample. 
The stars were then manually spectral typed with the HAM- 
MER pipeline (Covey etal. 2007); the error in the process 
is expected to be smaller than one sub-type, as discussed by 
lWestetalJ(l20TTl) . 

The radial velocities {Vr) of the stars were measured by 
cross-co rrelating the spectra using IDL routine XCORL.PRO 
dMohan tv & Basri 2003; West & Basri 2009 ) wifli the appro - 
priate spectral type templates from iBochanski et al.l ( I2007bl) . 
which are in the heliocentric rest frame. The cross-correlation 
was performed in the wavelength range of 6600-7550 A. 
Since templates are only available for MO-LO dwarfs, MO 
templates were used for K7 dwarfs. However, as each spectral 
type spans a range in mass and temperature, the spectral fea- 
tures in stars of the same subtype can differ significantly. This 
is probably the largest source of error in measuring Vr with 
the cross-correlation techniques. However, cross-correlation 
with templates typically yields Vr with < 10 km preci- 
sion for SDSS spectra (iBochanski et al.ll2007a) Other sources 
of error include difference in resolution between the template 
and object spectra and the accuracy of wavelength calibration, 
which was < 0.04 A/pixel rms (1.4 km s ^) for all but a few 
of the objects. 

To assess the errors in our radial velocities, we 
cross-correlated the observed radial velocity standards 

' As flux-calibrations were not taken for the first run, we used a subset of 
radial velocity standards, which had absolute flux measurements as part of 
the Palomar-Michigan State Survey (Reid et al. 1995), as flux standards. 

IRAF is distributed by the National Optical Astronomy Observatory, 
which is operated by the Association of Universities for Reseai'ch in Astron- 
omy (AURA) under cooperative agreement with the National Science Foun- 
dation. 



dPe lfosse et alj [l998h with appropriate templates from 
IBochanski et al.l mom . Figure [T] s hows the com p arison 
between our measured values and the iDelfosse et al.l ( Il998h 
values, which were measured from high-resolution spectra. 
Apart from three outliers from the 2010 seasons, our values 
compare well with the Delfosse et al. (1998). The median ab- 
solute deviation (m.a.d.) of the difference was 5.7 km s^^; we 
adopt 6 km s^^ as the typical error in our measurement of Vr- 
The spectra were then corrected for the measured radial ve- 
locities to be in the heliocentric rest frame and fed back into 
the HAMMER pipeline to measure the equivalent width of 
Ha; the molecular band strengths of CaH2, CaH3, and TiO 
features; and the S/N of the spectra, which was measured in 
the region spanning 6500-6550 A. 

3. RESULTS 

Figure |2] shows the spectral type and r — z color distribu- 
tions of the primary and the secondary components of the 113 
SLoWPoKES pairs that were observed (TableO; the number 
of pairs in each bin is also shown. The observed sample, ex- 
cluding pairs that were rejected for low S/N or other reasons, 
spans the K7-M4 spectral types (r — z = 0.66-2.52) for the 
primary and K7-M5 (r — z = 0.77-3.08) for the secondary. 
Even though our observed sample was limited to r ^ 17 and, 
thus, a dearth of late-type M dwarfs was to be expected, there 
are nonetheless 1 1 pairs with at least one component later than 
M4 and only two pairs with both components later than M4. 



3.1. Spectroscopic Binaries 

As discussed in Section |2] the radial velocities were mea- 
sured by cross-correlating the progr am spectra with the appro- 
priate SDSS template spectra (Bo chanski et alj|2007bh . The 
cross-correlation function (CCF) is used to determine the best 
match between the templates and program spectra. Pres- 
ence of multiple turning points in the CCF as well as un- 
usual broadening of the spectr um can indicate the p resence 
of a spectroscopic binary (SB; lMatiievic et al]|2010l) . While 
such a detection is unambiguous only in high-resolution spec- 
tra, we found possible evidence of SBs in our low-resolution 
spectra. Alternatively, the wide CCF could correspond to fast 
rotators, although our v sin i resolution of 35 km s^^ means 
they would have to be rotating at very high speeds. Figure [3] 
shows the CCF for the ten SB candidates (Table O. AH ten 
candidates have a relatively high S/N ratio, so the CCF is not 
a product of noisy spectral features. For comparison, the CCF 
for the radial velocity standards, which are presumably single 
stars, of the corresponding spectral types are shown in red, 
dashed lines. High-resolution spectra are required to confirm 
these SBs. 

Previous studies have found that components of wide bi- 
naries are more likely to have a companion as compared to 
single field stars. This enhanced binarity has been ascribed to 
the ease of transfer of angular momentum that facilitates the 
formation of close pa irs jTokovinin 1997; Bate et al. 2003; 
iBurgasser et alj 120051 ; IConnellev et al., ,2009.) and/or the sta- 
biUty of wide pairs in the field dLaw et all 2010 ). Among very 
low-mass wide binaries, the frequency of tight companions 
is (50±11)% tFahertv et al.. .2010.) . In a sample of nearby 



4 



Dhital et al. 



Table 2 

Properties of observed SLoWPoKES binaries 



ID 


Right Ascension 


Declination 


distance 


Ma 














A 


B 


A 


B 


A 


B 


A 


B 


A 


B 


A 


B 


SLW J2000 


(hh mm ss.s) 


(dd mm ss.ss) 


(pc) 


(mas yr 




(mas yr 




(km s 




1512+20 


15 12 22.52 


15 12 25.41 


+20 28 20.6 


+20 28 12.3 


63 


58 


-49 


-48 


4 


7 


-84.8 


-77.8 


0831+36 


08 31 23.12 


08 31 23.16 


+36 54 41.8 


+36 54 17.2 


70 


83 


49 


46 


-54 


-57 


-30.8 


-25.9 


0741+19 


07 41 55.34 


07 41 57.06 


+19 55 45.8 


+19 55 33.3 


66 


78 


-36 


-35 


-27 


-27 


31.2 


47.8 


0957+37 


09 57 57.18 


09 57 55.63 


+37 56 02.4 


+37 55 52.8 


87 


75 


-23 


-23 


-60 


-59 


-23.6 


-31.2 


1120+20 


11 20 03.38 


11 20 05.26 


+20 46 53.2 


+20 46 54.9 


96 


101 


-37 


-41 


-2 





-7.6 


-40.6 


0858+09 


08 58 57.80 


08 58 54.73 


+09 36 59.1 


+09 37 23.7 


65 


63 


-111 - 


107 


6 


6 


23.5 


32.7 


1527+49 


15 27 52.04 


15 27 50.57 


+49 08 54.2 


+49 09 47.4 


70 


65 


-60 


-63 


50 


53 


-89.7 


-90.2 


0734+28 


07 34 50.75 


07 34 49.23 


+28 17 39.7 


+28 18 15.8 


62 


75 


-27 


-29 


-28 


-29 


-2.3 


-36.7 


1318+47 


13 18 15.49 


13 18 15.00 


+47 30 29.4 


+47 31 33.7 


47 


47 


-104 - 


103 


33 


35 


-60.1 


-71.6 


1508+06 


15 08 44.07 


15 08 43.72 


+06 46 25.9 


+06 46 35.5 


110 


108 


-42 


-44 








-80.4 


-51.3 



Table 2 



ID 


r 




r 


— z 


Spectral Type 


EW (Ha) 


CaH2 


CaH3 


Ti05 




A 


B 


A 


B 


A 


B 


A 


B 


A 


B 


A 


B 


A 


B 


SLW J2000 


(mag) 










(A) 














1512+20 


15.01 


16.87 


2.02 


2.56 


M3 


M4 


-0.18 


0.00 


0.45 


0.36 


0.72 


0.63 


0.44 


0.34 


0831+36 


16.13 


16.51 


2.12 


2.36 


M3 


M4 


2.36 


7.93 


0.42 


0.38 


0.69 


0.62 


0.43 


0.37 


0741+19 


16.29 


16.67 


2.22 


2.46 


M3 


M4 


3.06 


7.18 


0.43 


0.41 


0.69 


0.67 


0.42 


0.38 


0957+37 


17.28 


17.34 


2.45 


2.57 


M3 


M4 


4.11 


7.16 


0.39 


0.39 


0.65 


0.64 


0.38 


0.34 


1 120+20 


16.62 


17.09 


2.13 


2.32 


M3 


M3 


7.16 


0.99 


0.46 


0.44 


0.72 


0.70 


0.45 


0.45 


0858+09 


15.58 


15.76 


2.12 


2.20 


M3 


M3 


5.12 


2.57 


0.42 


0.42 


0.65 


0.67 


0.39 


0.43 


1527+49 


16.00 


17.01 


2.19 


2.57 


M3 


M4 


4.83 


1.08 


0.32 


0.41 


0.59 


0.68 


0.31 


0.41 


0734+28 


14.40 


14.78 


1.65 


1.90 


M2 


M2 


4.64 


0.00 


0.46 


0.53 


0.69 


0.76 


0.49 


0.55 


1318+47 


15.19 


17.40 


2.20 


2.95 


M2 


M5 


1.95 


3.94 


0.41 


0.33 


0.68 


0.62 


0.42 


0.27 


1508+06 


15.66 


17.04 


1.78 


2.23 


M2 


M2 


3.53 


0.64 


0.42 


0.49 


0.67 


0.75 


0.43 


0.52 



Note. — The first 10 pairs are listed here; the full version of the table is available online. 





M5 








M4 






(J 


M3 






o 

a, 


M2 




1 




Ml 




5 




MO 




7 




K7 




3 



1 5 2 



^^^^^^ 



8 6 



3.0 














3 




2.6 












1 


1 


1 1 














1 


3 


2 


2.2 






1 




1 




6 5 


2 














7 


9 10 






1.8 








2 


2 




6^ 










1 




6 


1 








1.4 




1 


3 


6 


1 












- 1 




5 












1.0 




1 
















- 1 












1 


1 1- 


0.6 








1-5 


5-10 




10-15 15-10 



K7 MO Ml M2 M3 M4 M5 



0.6 1.0 1.4 1.8 2.2 2.6 

(r-z)p„ 



3.0 



Figure 2. The distribution of measured spectral types and r — z colors for both the primary and secondary components of the 113 SLoWPoKES 
pairs observed in our KPNO run. The spectral types are shown in 1 subtype bins while the colors are in 0.2 mag bins. The background colors 
indicated the density in each bin, with the individual numbers printed as well. By definition, the primary always has a earlier spectral type and 
a bluer color than the secondary. 

hierarchical, we were probi ng a diff erent kind of hierarchical 
systems than those found bv lLaw et al. 2010. Here we probed 
the extremely close pairs (spectroscopic binaries) whereas 



SLoWPoKES pairs, iLaw et al.l (|2010^ found that the bias- 
corrected higher-order multiplicity was While only 
10 of 113 pairs (8.8%) in this study have been identified as 
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Table 3 

Candidate Spectroscopic Binaries 



Name 


Spectral Type 


Note 


SLW 1211+58a 


MO 




SLW 1320+43b 


Ml 


wide CCF 


SLW 1417+13b 


Ml 


wide CCF 


SLW 1441-02a 


K7 




SLW 1340+27a 


Ml 


wide lines; wide CCF 


SLW 1429+42a 


MO 


double peak 


SLW 1508+06b 


M2 


wide CCF 


SLW 0848+26b 


Ml 


wide CCF 


SLW 0934+15a 


M2 


line splitting; double peak 


SLW 0734+28a 


M 


line splitting, wide CCF 



iLaw et ai] |2010 probed systems with separations l arger than 
8-10 AU. In fact, our results are consistent with the lLaw et alj 
120101 findings but limited to the extremely close pairs. 



3.2. Fidelity of SLoWPoKES Pairs 

The observed pairs were identified in iPaper II based on a 
matching of their position, distance, and proper motions. The 
third velocity component, Vr, can be used to test the fidelity of 
the observed pairs and, by extrapolation, of the SLoWPoKES 
catalog. 

Figure |4] shows the radial velocities of the primary compo- 
nent against that of the secondary in the left panel and the 
distribution of their differences in the right panel. The iden- 
tified candidate SBs are shown as concentric circles. Exclud- 
ing the ten candidate SBs, 90 of the remaining 103 pairs (i.e., 
87.4%) have AVr within 3 <t of the mean; the 3 a region is 
shown in gray in the left panel. Overall, the AVi distribu- 
tion is well-fit by a Gaussian with ji = —0.97 ± 0.80 km s~^ 
and a ~ 12.04 ± 0.80 km s^^, shown as the dashed line 
though there are more outliers than expected (see below). 
As a w 1.4826 m.a.d for large normally distributed popula- 
tions, (TAVr ( = X ^v-r — "/2 X 1.4826X 5.7 km s^^ = 
11.95 km s^^) is in excellent agreement with the m.a.d. we 
measured for our radial velocity standards. 

We checked that the AVr distribution of our binaries is in- 
deed distinct from ph ysically unassociated stars. We used the 
Galactic model from iPaper II that gives the expected 3-d ve- 
locity distribution for any position in the Galaxy or, if de- 
sired, a randomly chosen velocity from that distribution. We 
compared with a sample that has been selected in a similar 
manner to the pairs in our observed sample. Hence, at the 
Galactic positions of each of the observed pairs, we gener- 
ated pairs of 3-d velocities until a pair with matching proper 
motions was found. The proper motion matching criteria was 
the same as that in Paper 1. For statistical robustness, we con- 
ducted 10*^ realizations of this simulation; the normalized his- 
togram of the resultant distribution is shown in dotted lines 
in the right panel of Figure |4] For a quantitative assess- 
ment of the difference between the simulated and observed 
AV r distributions, w e performed the Kolmogorov-Smirnov 
test (" Press et al.ll 19921) and found a 0.93% probability that the 
two were drawn from the same parent population. We con- 
clude the AVr distribution of our observed binaries is much 
narrower than the scatter expected of two unassociated stars. 

Thirteen (12.6%) of the pairs have Vr that disagree at >3(t. 
The SLoWPoKES catalog only contains pairs with probabil- 
ity of chance alignment, Pf , less than 5%, meaning fewer than 



five of 103 pairs were expected to be fal se positives. In fact, 
Pf was tabulated for each pair in IPaper 3 ; the cumulative sum 
of chance alignments was only 0.3%, implying <1 pair was 
expected to be false positives. As shown in the left panel of 
Figure |5] there is no trend in AVr as a function of the Pf . 
There are discrepant pairs at all values of Pf. However, as 
can be seen in the right panel of Figure |5] there is a signif- 
icant trend of AVr with the S/N in the spectra. The pairs 
with the largest AVr values are at low S/N while there are 
no discrepant pairs at high S/N. This suggests that the cross- 
correlation process and the measurements of Vr might have 
been adversely affected by the noise, yielding noisier radial 
velocities. In addition, given the large hi gher-order multi plic- 
ity fraction seen in SLoWPoKES (45 % ; Law etani20 1 Oh . the 
presence of more SBs in our observed sample cannot be ruled 
out by our low-resolution spectra. 

In summary, the vast majority of the sample pairs show 
agreement in their radial velocities, as expected for physi- 
cal binaries. Pairs with discrepant Vr^ have spectra with low 
observed S/N or are (candidate) hierarchical systems with a 
spectroscopic binary. 

3.3. Metallicity Index Calibration 

Figure |6] shows the observed SLoWPoKES pairs, with 
the components of each pair connected with a solid line, 
in the CaH-TiO space w ith the dotted lines showing the 
iso-C lines from ILRS07I for C, = 0.1-1.5 in steps of 
0.1. The solid lines delineate the boundary between the 
dM/sdM/esdM/usdM classes (C 0.825, 0.500, and 0.200, 
respectively); the single-star spectral standards for the sdM 
(diamonds), esdM (triangles), a nd usdM (squares) classes are 
also shown (iLepine et al.l 2007b . For clarity, only pairs whose 
error bars, in both CaH and Ti05 of both pairs, are smaller 
than the median error are plotted. Most of the observed pairs 
are dMs, i.e., part of the thin disk with roughly solar metallic- 
ity. This is not surprising for a bright sample located within 
^•200-300 pc of the Sun as t he local neighborhood is largely 
dominated by thin disk stars (iBochanski et al.ll2010l) . 

Most of the binary pairs lie parallel to the iso-^ lines, within 
the error bars. The inset in Figure |6] shows the distribution of 
AC,, which is centered around zero but has substantial scat- 
ter Quantitatively, the median and median absolute deviation 
(m.a.d.) are -0.005 and 0.066, respectively, indicating that for 
the vast majority of the observed pairs, C is a correctly infers 
a common metallicity for the two stars in each pair 

To test whether the observed A^ distribution was merely 
the intrinsic scatter in A(, we randomly selected 113 pairs of 
disk stars (C, = 0.825-1.100) in the similar brightness range 
as the SLoWPoKES sample (r < 20) and with high proper 
motions (/n > 40 mas y r~^) from the SDSS spectroscopic 
catalog dWest et al.ll201 Ih and calculated the A^ distribution. 
There were 8030 stars in the DR7 spectroscopic sample that 
met these criteria; 113 pairs were randomly selected from this 
sample and their distribution calculated. We performed 
this simulation 10^ times; the resultant distribution is plot- 
ted in red in the inset of Figure |6] The simulated A^ dis- 
tribution is less centrally peaked, and much broader, com- 
pared to our observed distribution. The Kolmogorov-Smirnov 
test (Press et al. 1992) finds a probability of 3.3x10"^^ that 
the two distributions were drawn from the same parent pop- 
ulation. We conclude that the metallicity of components of 
SLoWPoKES binaries observed in this program are more sim- 
ilar to each other than that of two randomly paired thin disk 
stars. As components of a binary system are expected to have 




Figure 3. Broad cross-correlation functions (solid lines) and/or split spectral features indicate the presence of a tight companion in the ten of the 
observed systems. For reference, the auto-correlation function of our radial velocity standards are also shown (dashed lines). High-resolution 
spectra are needed to confirm the spectroscopic binaries. All spectra were corrected to the heliocentric rest frame. 
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Figure 4. Left: The radial velocity of the primary and the secondary components: the pairs with AVr > 3ct are shown as triangles while the 
candidate spectroscopic binaries, shown as red concentric circles, have larger AVr. Right: The distribution of difference in radial velocities 
between components of SLoWPoKES pairs observed in this program, with the Gaussian fit shown in dashed lines. 85% of the sample has AVr 
< 3(T. All pairs with AVr>50 km s^ have relatively low S/N.. The expected intrinsic scatter in AVr for two una ssociated stars in the same 
3-d position in the Galaxy and have matching proper motions, as calculated using the Galactic model from iPaper H is shown with dotted lines; 
it is a much larger dispersion compared to our observed sample of CPM pairs. 



formed of the same material, this further serves to confirm the 
physical association of the pairs. Furthermore, it strengthens 
the argument that is a reliable proxy for relative metallicity. 

However, Figure |6] also demonstrates some deficiencies 
in the definition of C,. First, /S.C, is more than three m.a.d. 
away from zero for ^18% of the pairs, versus the ^5% ex- 
pected for a normal distribution. They are discrepant es- 
pecially at large values of Ti05 and CaH2+CaH3, i.e., for 
higher-mass M dwarfs, perhaps suggesting a break in the 



C(Ti05, CaH2 + CaH3) relation. Large errors in this regime 
further complicate the issue, as the discrepant C, values could 
result from the difficulty of measuring the shallower Ti05, 
CaH2, and CaH3 band heads in late-K and early-M dwarfs. It 
is also evident how the iso-^ contours converge at the higher 
masses as pointed out out by LRS07. On the other hand, the 
discrepancy persists for the higher-mass pairs with smaller er- 
ror bars, as can be clearly seen in Figure |6] Second, and per- 
haps more importantly, our measured C, values increase and 
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Figure 5. Left: The probability of chance ahgnment (Pf) calculated using the Galactic model by 'Paper T vs. the difference in radial velocity 
between components for 113 SLoWPoKES CPM pairs observed in this program. Right: Systems with a lower S/N (red) have a higher dispersion 
in AVr, suggesting the observed discrepancy might be due to the lower S/N. Candidate SBs are shown as red, concentric circles while pairs 
with AVr > 3a" are shown as triangles. 
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Figure 6. Distribution of the observed M dwarf binaries, with components of a system connected by solid lines, in the CaH/TiO space confirms 
the iso-( contours with ( — 0.0, 0.1, 0.2, 1.5 (dotted lines), albeit with some deviation at the highest values (towards top right). The 
distribution of A( (blue histogram; inset) peaks at ~0 and is significantly different from a distribution for two randomly paired dMs (red 
solid line). However, the scatter in A(" is large; the bar at the middle-top of the inset shows the scatter in the error of A(J (itac)- The solid 
contours demarcate the boundary between the dM/ sdM/esd M/usdM classes while the K7-M5 spectral standards for the sdM (diamonds), esdM 
(triangles), and usdM (squares) classes are shown l lLRSOTh . Early-type M stars are at the upper right, late-types at the lower left. 

120081) . Yet, the ^-index yields the opposite. This result neces- 
sitates a redefinition of (. With a spectroscopic sample of 



become super-solar (i.e., C > 1) for the higher mass stars. 
This is inconsistent with the expectation: given the apparent 
magnitude constraints (r ^ 15-17), the higher mass stars in 
this sample can be expected to be farther away and, hence, at 
larger Galactic heights given most of the SPSS sight lines are 
at at high Galactic latitudes ("I vezic et al .''20b8). Stars at high 
latitudes are, on average, older; and consequently, i f anything, 
they might be expected to have lower metallicities (IWest et al.l 



113 visually resolved binaries, we are in an unique position to 
modify the definition of 

Given the lack of subdwarf pairs in our sample, any recal- 
ibration of ( would be systematically biased to high metal- 
licity. Hence, we conducted a search for companions around 
the subdwarfs (C < 0.825) in the SDSS DR7 spectroscopic 
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Table 4 

Coefficients, a jv , for Eq J2j 



Coefficients 


LRS07 


This study 


ao 


- 0.050 


- 0.047 


ai 


- 0.118 


- 0.127 


a2 


+ 0.670 


+ 0.694 


as 


- 0.164 


- 0.183 


34 




- 0.005 



catalog jWest et al.ll2011l) with extant SDSS spectra. We have 
identified a sample of ten pairs with low values of (; and they 
are shown as open circles in Figure |6] The full sample that 
was identified in the search will be presented in a future pa- 
per (Dhital et. al., in prep.). One of the added pairs is at the 
sdM/esdM boundary while the other nine are at the dM/sdM 
boundary. While small this sample provides an invaluable 
constraint in the low-metallicity regime. 



|LRS03definedCas: 

C= (1) 

^ 1 - [Ti05]zg ' ^ ' 

where [Ti05]zQ is a third-order polynomial of 
(CaH2H-CaH3): 

[Ti05]zo = H aJv(CaH2 + CaH3)^, (2) 

N 

and where the coefficients, aN, are tabulated in Tableland 
were obtained as a single fit to the Ti05 index as a function of 
CaH2H-CaH3 index for kinematically-selected sample of thin 
disk stars. 

We can recalibrate C by varying the functional form of 
[Ti05]zQ in Eq. Q such that the scatter in the A( distri- 
bution is minimized and distributed around zero. As noted 
earlier, the LRS07 definition, to first-order, is a robust mea- 
sure of relative metallicity; and a recalibration need only be 
a perturbation about that definition. Moreover, as the defini- 
tion was based on the distribution of (Ti05, CaH2H-CaH3) of 
disk stars, it is a good starting point for the recalibration. So 
we on ly chose to explore the coefficient values within ±0.03 
of the Lep ine et all (12007!) values in steps of A = 0.001. We 
have introduced a fourth-order term Eq. (|2j, with an initial 
guess of zero, based on the observed deviation of higher-mass 
pairs from the iso-( lines. The best fit values for the coef- 
ficients were found by minimizing x^, where = was 
assumed to be the model. All dM/sdM pairs, except for the 
ones with the large error bars, were considered for the fit. 

Figure|7]shows the new iso-( contours, with the coefficients 
tabulated in Table|4] The contours look significantly different 
despite small changes in the coefficients, reflecting the very 
sensitive dependence of ( on its independent variables. The 
differences can be summarized as: 

(i) The scatter in the A^ is smaller, with the m.a.d. de- 
creasing from 0.060 to 0.044. There are fewer outliers 
as well, with the values converging towards A( — 0. 
As the median error in is bigger than the scatter in 
A(, decreasing the scatter further is not possible un- 
less higher S/N data are obtained. The A^ distribution 
is much more centrally peaked compared to randomly 
paired field stars. 



(ii) With only a few pairs at the higher-masses (large Ti05), 
C yields more consistent values between components of 
the higher-mass dMs than in the original iLRS07 defi- 
nition. This congruence is important if we are to have 
the same metallicity proxy for all low-mass stars and is 
likely to improve as more pairs are added to that locus. 
It is especially reassuring to get the same value of ( for 
the components of pairs with large differences in mass, 
CaH2H-CaH3, and Ti05; the previous definition such 
pairs were especially discrepant as the primary was as- 
signed a artificially super-solar (. With the new defi- 
nition, the higher-mass dMs instead have slightly sub- 
solar (, which is what is expected for stars at larger 
Galactic heights. 

(iii) The new iso-( contours are less crowded at the higher- 
mass regime, allowing for a more robust determina- 
tion of metallicity. In fact, the iso-( do not converge 
before (Ti05, CaH2H-CaH3) = (1,2), hence, expand- 
ing the regime for sdMs/esdMs/usdMs. The metallic- 
ity classes also become more sensitive to Ti05 rela- 
tive to CaH2H-CaH3. However, with the new contours 
for C the previously defined standards for the metallic- 
ity classes are assessed to be more metal-poor and no 
longer are in the same class. While this argues for a 
definition of new standards, we advise against such a 
revision until there are more subdwarf binaries to more 
robustly calibrate the contours in that regime. 

In general, the new definition of ( better fits the observed 
sample of visual binaries as well as resolving outstanding is- 
sues at the high-mass end. However, due to a paucity of sub- 
dwarf pairs, it leaves the low-metallicity regime rather uncon- 
strained. ( looks to be a good proxy for metallicity, and future 
observations of subdwarf pairs should calibrate it for all low- 
mass stars. Studies that are using subdwarf binaries (Dhital et 
al., in prep.) and F/G-K/M binaries (Bochanski et al. in prep.; 
Lepine et al. in prep.) are akeady underway and are part of a 
larger effort to measure the absolute metallicity of low-mass 
stars. 

3.4. Ha Activity 

As M dwarf photospheres are too cold to excite elec- 
trons into the n = 2 level, any observed Ha feature, ei- 
ther in absorption or in emission, is chromospheric in ori- 
gin. The weakest chromospheres will exhibit no Ha; as the 
activity levels increase. Ha will be observed in absorption 
with the line filling in and eve ntually going into emission 
for the most active M dwarfs dsFauffer & Hartmann| 119861 : 
ICram & Giariip^ [19871: IWalkowicz & Hawlevll2009h . As 
Ha is in absorption for both inactive and moderately active M 
dwarfs. Ha emission has traditionally been used as the tracer 
of chromospheric a ctivity and is bias ed towards the most ac- 
tive M dwarfs (e.g.. lWest et al.ll201 Ih . 

In our sample, 11 of the 113 pairs showed clear evidence of 
Ha emission in both components while three pairs had only 
one component with Ha in emission. In addition, 33 pairs 
showed Ha absorption in both components and 22 pairs in one 
component; the remaining 47 were classified as inactive. The 
large fraction (39%) of stars with Ha absorption is in accor- 
dance with the nearby M dwarfs in the Palomar/MSU Nearby 
Star Survey Spectroscopic Survey (iGizis et al.l 120021) . Over- 
all, for 91 of the 113 (81%) pairs in our sample, both compo- 
nents of a pair showed similar levels of activity — in emission. 
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Figure 7. Same as Figur ed but with t he red efined (, which was done by adding a fourth-order term in Eq. l|2} and refitting the coefficients 
by perturbing about the iLepine et alj ( 120071) values. The new definition of minimizes the scatter in the A( distribution and yields lower 
metallicity for stars at larger Galactic heights, which are expected to be older and have a lower metallicity. The Ai^ distribution is now also 
much more centrally peaked compared to random pairings of unassociated stars (compare to inset of Figure|6] 



absorption, or the lack or activity. For the pairs with both 
components having Ha emission, we converted the equiv- 
alent width in Ha to Lna/Lboi — the ratio of Ha luminos- 
ity to the bolometric luminosity that is independe nt of spec- 
tral type — following iWalkowicz et al.l (120041) and IWest et all 
(12004) : all pairs had comparable levels of Lna/Lboi within 
the error bars except for the two where the active primary was 
identified as a candidate SB and had an inactive secondary. 
The tidal forces due to the tight comp anion has presu mably 
enhanced th e activity of the primary (IShkolnik et al.l [20101: 
ISilvestri et al. 2006). 

As Ha activity depends strongly on mass (IWest e t al. 2008), 
one way to compare the intrinsic variability in activity lev- 
els is by only looking at pairs with components of similar 
masses. All sixteen pairs with Ha in emission or absorption 
in both components and with similar colors (A(r — z) < 0.2) 
had Ha equivalent widths within 130% of each other Com- 
pared to the 200-300% difference in H a activity exhibited 
by M dwarfs over time (iBell et al.ll201 ll) . the components of 
the similar-mass binary pairs in our sample exhibit a much 
smaller difference in Ha activity. This is consistent with the 
expectation that stars of similar ages and masses have com- 
parable activity levels, presumably because they experience 
similar spin-down rates. 

4. CONCLUSIONS 

We have carried out a spectroscopic follow-up study of 11 3 
ultra-wide, low-mass CPM binary systems from the SLoW- 
PoKES catalog (.Paper L) using the GOLDCAM spectrograph 



on the KPNO 2.1 m telescope. We measured the radial veloc- 
ities of each component by cross-correlating them with ap- 
propriate standards and used them to assess the fidelity of 
pairs in the SLoWPoKES catalog. 95 of the 113 (84%) of 
the pairs have the same radial velocity within 3 a. At least 
five of the pairs with discrepant radial velocities are candi- 
date SBs, which would explain the difference. There may 
be additional spectrosc opic comp anions undetected in our 
low-resolution spectra. iLaw et ail (2010) found that 45% of 
the SLoWPoKES systems are either hierarchical triples or 
quadruples. Either high-resolution spectroscopy or imaging 
would be needed to identify the close companions and to fur- 
ther quantify the incidence of higher-order systems in wide 
binaries. 

We examined the Ha activity in our observed sample. The 
components of binary pairs exhibited overwhelmingly com- 
parable levels of Ha activity. Moreover, the AHa of the pairs 
with similar r — z colors and two active components, while 
large, was several times smaller than the variation seen in sin- 
gle M dwarfs. Our results coiToborate that low-mass stars of 
the same mass should spin-down at the similar rates over time. 
However, larger samples of active are needed to confirm this 
finding and to constrain the rate of this spin-down. 

We tested the LRS07 ^-index and found that, to first-order, 
it is a robust measure of relative metallicity. The value of ( for 
the two components in each binary system match within the 
error bars for most pairs, indicating a common metallicity as 
expected. However, we find a systematic bias for the higher- 
mass M dwarfs such that ( overestimates the metallicity. As- 
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suming all of the pairs are physically associated systems and 
have the same metallicity, we have redefined While the 
shift is small, it better represents iso-metallicity lines in the 
high-metallicity regime and represents an incremental step 
towards defining an absolute metallicity scale for low-mass 
dwarfs. Planned further observations should extend the cali- 
bration of C as well as map it to an absolute metallicity scale 
in the near future. 
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